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Abstract. 


Large water movements do exist in frozen soils. They are induced by a depression 
on the frost line (isotherm 0°C) or in the frozen zone of the soil. This depres=- 


sion is the cryogenic suction. Those movements of water are the main causes for 
both frost heave and damages to structures built on or inside a frozen soil. 


We describe a model computing both temperature and hydraulic potential. Practical 
situations are investigated through numerical simulations. 


INTRODUCTION 


Une des principales difficultés du dimensionnement des 
structures (routes, gazoducs, batiments,..) dans les 
régions soumises réguliérement 4 des températures né=- 
gatives est l'évaluation de la gélivité du sol qui 
doit supporter ces structures. La gélévité est une no- 
tion globale qui permet de classer les sols suivant 
leur comportement lorsqu'ils gélent [1], [2], [3]. Les 
sols humides gélifs ont leurs caractéristiques consi- 
dérablement modifiées par le gel ; ils subissent un 
gonflement qui peut étre notable (parfois 9.5m par an 
([4]) et donc trés nuisible aux structures qui reposent 
sur ces sols. Ce gonflement est di a une forte aspira- 
tion d'eau vers le front de gel et les zones gelées : 
c'est la succion cryogénique. Elle provoque d'impor- 
tants mouvements d'eau vers le front de gel et une ac- 
cumulation de glace dans les zones gelées, ce qui peut 
donner naissance a de fortes contraintes et déforma- 
tions dans le sol. 


Nous présentons une modélisation des comportements 
thermique et hydraulique des sols gelés qui décrit la 
succion cryogénique. Le modéle, basé sur des notions 
macroscopiques usuelles de la mécanique des milieux 
continus, permet de calculer 1'évolution de la tempé- 
rature et de la charge hydraulique (et donc les mouve- 
ments de l'eau dans le sol) par une méthode d'éléments 
finis. 


MOUVEMENTS D'EAU DANS LA CONGELATION DES SOLS 


La quantité d'eau contenue dans un sol et les mouve- 
ments de cette eau sont des aspects primordiaux de 
l'étude de la congélation d'un sol. En particulier 
l'existence d'eau liquide dans les zones gelées (méme 
a des températures de plusieurs degrés C inférieures 
4a °C) a une grande importance dans le comportement 
thermique et hydraulique du sol [5], [6]. On explique 
généralement la présence d'eau liquide dans les sols 
gelés par des arguments de chimie et physique micros- 
copique (par exemple des phénoménes de capillarité et 
d'adsorption) [6]. 


Nous supposons dans notre modéle, comme cela est sou- 
vent constaté expérimentalement, que pour un sol donné 
la proportion volumique d'eau liquide dans le mélange 
eau liquide-glace ne dépend que de la température @ 
dans le sol (2), [5], [6]. Cette proportion volumique 
v(8) est mesurée expérimentalement et une fonction 
v(8) est représentée sur la figure l. 


V(®) 


1 


discontinuite de V 
—> 


= S3C 0 5°C 


Figure l 


Concentration volumique d'eau liquide 
dans le mélange eau liquide-glace. 


On constate également que, bien que la perméabilité 
du sol gelé soit trés faible, l'eau liquide se déplace 
dans le sol gelé [4]. Nous supposerons que la perméa- 
bilité du sol est une fonction de la température [4]. 


MODELISATION DES PROBLEMES THERMIQUE ET HYDRAULIQUE 
Nous nous placons dans le cadre de la mécanique deg 
milieux continus pour écrire les équations de conser- 


vation et les lois de comportement [7]. 


Nous supposons que le sol est un milieu poreux saturé, 
il y a donc trois constituants : le sol, l'eau liquide 
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et la glace. Les inconmes du modéle sont la tempéra- 
ture et la vitesse de l'eau en tous les points du sol 
et & chaque instant. Nous ne nous intéressons pas ici 
aux mouvements du sol et nous sommes conduits & suppo- 


ser connue en tout point et 4 tout instant, la porosité 


€ du milieu, qui représente la proportion du volume 
d'eau et de glace dans le volume total. Un modéle plus 
complet incluant la prise en compte et le calcul des 
mouvements du sol est actuellement en développement 


{13]. 


Nous écrivons l'équation de conservation de la masse 
d'eau (liquide + glace) et 1'équation de conservation 
de 1'énergie. 


L'équation de conservation de la masse d'eau s‘écrit : 


(1) 


at o+ eh =0Q, 


~ 


ou : 


Yk eve + e(l-v)o , Pe. et ee étant les masses volu- 


miques respectives de l'eau et de la glace (kgem73) ; 


A me est la vitesse de l'eau (m*s~!). 
L'équation de conservation de l’énergie s‘écrit en né- 
gligeant les mouvements du sol : 


(2) 


om div(e U.) +divq#tr, 


ae =.CO+p eve, 
aC mec cect cable eg SIC,” 
- e wp ev 6+9 eEvk . 
e ee e 
. 6 est la température (°K), 


. & est la chaleur latente massique de l'eau 


(J3.kg73), 


~ 1Gey Pe Ea sont les capacités calorifiques massi- 


e 
ques de l'eau, de la glace et du sol 


teks" 2.K3"); 


» 0, est la masse volumigue du sol, 


. g est le vecteur courant de chaleur (J.m~?.K7?), 


2 1 


- © est le taux volumique de chaleur_regu J.m~“.s™", 
La fonction v (figure 1) donne 4 l'énergie interne du 

milieu e la forme représentée sur la figure 2. Sur les 
fronts de gel (@=O0°C) il faut donc considérer 1'équa- 


tion (2) comme une équation aux discontinuités. 


En portant l'équation (1) dans 1'équation (2) et en 
négligeant un terme de convection d'énergie calorifi- 
uve transportée par le mouvement de l'eau 
fo.ew,grad(C,9)), on obtient pour la température 6, 


une équation ot la vitesse ee n'intervient plus 


5 Cc 2Seren arnt 


at g (3) 


+divq®= r. 


Pour les lois de comportement, nous choisissons de 
maniére classique : 


Figure 2 
Energie interne du milieu fonction de la température 


— la loi de Fourier 


q= -\ grade, (4) 
ou 4 est la conductivité thermique du milieu 
(3. 872.05" Ks). Elle peut dépendre de la température. 


— la loi de Darcy 
m 


e 
ue ~ Zp eradh, 


(5) 


ou h est la charge hydraulique (h=p+ p82 ou p est 
la pression d'eau), m est la mobilité du milieu 
(m3.s.kg-!) relige & la perméabilité K par m= 7 


avec p= 1073 (viscosité dynamique de l'eau). 


L'équation (3) constitue un probléme 4 frontiére libre 
la fonction (l-v)e étant discontinue sur les fronts 
de gel a priori inconnus (puisqu'ils dépendent de la 
solution 6). 


Pour étudier ce probléme, nous intégrons l'‘équation (3) 
par rapport au temps et nous introduisons une nouvelle 
variable : l'indice de gel [8]. 


Pour simplifier, on suppose ici que C et A sont des 
constantes (voir [9] pour le cas général). En inté- 
grant (3) par rapport au temps, on obtient : 

EN t 
6) = | r]=R, 
o/ Yo 
On définit alors l'indice de gel, bien connu des ingé- 


nieurs pour apprécier dans certains cas la profondeur 
de gel (formule de Portnov [10)). 


cé- ae an se - div (ered r (6) 


t 
(x,t) = | 6(x,t)dt, 
0 
et l'équation (6) devient 
C 31 ( 1 aI 
" st oe rag | 
ou Al We iiys. se est le PETS As de 1. 


Ayant déterminé I par l'équation (7), la température 6 
est donnée par 


(8) 
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RESULTATS NUMERIQUES 


Plutot que d'exposer en détail la résolution par un 
schéma en temps semi-explicite ou implicite et par 
éléments finis de l'équation (7), nous préférons pré- 
senter des exemples de calcul en renvoyant a [9], [11], 
pour une étude compléte. Indiquons juste que connais- 
sant I on calcule 6 par (8), puis la charge hydrauli- 
que h grace aux équations (1) et (5). 


Le premier exemple que nous présentons illustre 1'in- 
fluence thermique locale d'un tube en aluminium enfoncé 
dans un sol refroidi par sa surface (figure 3). 


m 


Figure 3 


Tube d'aluminium dans un sol 


Bien sir, les conductivités ek du sol et de 
—. 140) - Nous sup- 
\s 

sons que le sol est uniformément saturé en eau €#= .2 

(bien stir ¢€=0O dans l'aluminium). Les données physi- 

ques relatives 4 cet exemple sont les suivantes 


(en SI) : 


l'aluminium sont trés différentes 


SUR 2000 dans le sol, 
ae. 2560 dans 1'aluminiun, 
Se 1000, 
7p =. 900, 
& 
» 2 = .334x108, 
- C_ = 4180, 
e 
- C_= 2090, 
g 
- C = 836, 
8 
- 4 = 1.42 dans le sol gelé, 
'. \ = 1.66 dans le sol gelé, 


» 4 = 200 dans 1l'aluminium. 


A l‘instant initial, on suppose que la température 
dans le sol est uniformément égale a 0.1°C. 


Les conditions limites sont les suivantes : 
— les cétés AD, CB et CD sont adiabatiques ; 


— la température &@ la surface du sol AB est refroi- 
die suivant la loi: 


6=0.1 t pour ¢%3600s, 


Sugg 
3600 
6= =2.9°C pour t > 3600s. 
Les figures 4 et 5 montrent l'état des températures 
dans le sol aux instants t= 300s et t=7200 s. On 
peut constater que la température dans 1‘aluminium 
tend rapidement a étre uniforme et égale a la tempé- 
rature de surface, et que le tube fait geler une par- 
tie du sol qui l‘entoure. La zone gelée est hachurée. 


01m 


Figure 4 
Isothermes aprés 5 minutes. Plus de la moitié de L'alu- 
minium est déja a température négative. La zone a tem- 
pérative est hachurée. 


Figure 5 


Isothermes aprés 2 heures. La température dans 
l'aluminium est constante. 
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Le deuxiéme exemple décrit -le comportement thermique 
et hydraulique d'un sol traversé par un gazoduc froid. 
Ce calcul a été effectué dans le cadre d'un contrat de 


recherche liant l'Université Carleton d'Ottawa, 1'E- A 
cole Nationale des Ponts et Chaussées (Paris), le Cen- Se VPS SS 
tre National de la Recherche Scientifique (Paris) et Ss A PISS 
le Laboratoire Central des Ponts et Chaussées (Paris). St SK : LSS 
Ce programme, parallélement au travail de modélisation SS TS >< >< 
zs ee . Ke Se 
et aux expériences en laboratoire, comporte une expé- ates Pi 1 >< 
rience en vraie grandeur concernant un gazoduc froid sts Seorseas 
enterré dans une station expérimentale (voir [1],[12]). aera Se SS 
L'expérience numérique que nous présentons utilise cer- aes Wisi s 
taines données physiques mesurées lors de cette expé- aS SS BR WW Seas 
rience (en SI) SSA PS SR TS 


[’ 


7 


S 

a 

Ne 25 
SSA 

man 

iS 

\ 


| WWW, 
eS ZINE NINN NI? 
- 4 =1.9 dans le sol non gelé, INANE IAAI 


= ans le so elé . N \ IN 
eA 202. 1 1 gelé, OR AANNNYN 
5 (@ ye ans le sol non gelé, VIN NAM /NV/N NY 
a ee ZZ PNP NS Ni 


» € = .39 dans le sol non gelé, 


- & = .334x108, 


+ COOE +00 nA 4 (S 22008001 OAR 
. C = 836, 
r 
. m © 107!) dans le sol non gelé, 
Fipure 6 


. m © 107!" dans le sol gelé. Maillage du sol autour d'un gazoduc 
La température initiale du sol est linéaire entre la 

base 8(y=0)=5°C et la surface du sol 6(y=2)#0.1°C. 

On calcule la température et la charge hydraulique 

toutes les deux heures. La figure 6 montre le maillage 

du sol entourant le gazoduc. 


La surface du sol reste constamment 4 la température HR 216000£ +02 
0.1°C. La base du sol (y=a) reste a la température 
5°C. La température 4 la surface du tube (cercle de la 
figure 6) est abaissée de la température initiale a 
-2°C. Les bords x=0O et x=2 sont adiabatiques. 


Pour la charge hydraulique, on suppose que la base 
(y=0) est alimentée en eau (h=0, la charge est défi- 
nie 4 une constante prés), alors que toutes les autres 
frontiéres sont imperméables (grad h.n=0). 


Les figures 7,8 et 9 présentent respectivement les 

isothermes, les profils de charge hydraulique suivant 
AA, BB et CC (figure 6) et les vitesses d'eau aprés ISOTHERNES 
16 heures. On peut constater, en accord avec l'‘expé- 


rience que l'eau monte de la base du sol vers le haut eee te 
du tube (figure 9) od il y a un front de gel (voir seg 
l'isotherme O°C sur la figure 7). La figure 8 montre Ase 


la dépression de l'eau sur le front de gel : elle se 
traduit par une forte variation de la charge d'eau : 
c'est la succion cryogénique. 


Les figures 10, 11 et 12 représentent les mémes quan- BOOOE= OORT 
tités (que 7, 8 et 9) aprés 48 heures. On peut noter Titoce o 
en particulier la progression du front de gel sur la BASE ip 
figure 10, et l'’aspiration de l'eau qui se trouve a 

la base sur la figure 12. 


Figure 7 


Enfin, la figure 13 montre les profils de température Isothermes dans le sol aprés-i68 


suivant AA, BB et CC a 48 heures. 
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Figure 8 
Profils de charge hydraulique aprés 16H 
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Figure 10 
Isothermes dans le sol aprés 48H 
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Figure 1l 
Profils de charge hydraulique aprés 48H. 
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Figure 12 
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CONCLUSION 


Notre modéle est basé sur des variables macroscopiques 
accessibles & l'expérience. Notre programme de calcul, 
disponible au Laboratoire Central des Ponts et Chaus- 
sées, peut @tre utilisé pour des problémes industriels 
d'ingénierie. I] permet de prévoir les comportements 
thermique et hydraulique des sols gelés 4 partir de 
données accessibles & 1'expérience. 
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